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ABSTRACT 

We assume that the relativistic sources moving along the dipolar magnetic 
field lines emit curvature radiation. The beamed emission occurs in the direction 
of tangents to the field lines, and to receive it, the sight line must align with the 
tangent within the beaming angle I/7, where 7 is the particle Lorentz factor. By 
solving the viewing geometry in an inclined and rotating dipole magnetic field, we 
show that, at any given pulse phase, observer tends to receive radiation only from 
the specific heights allowed by the geometry. We find outer conal components are 
emitted at higher altitudes compared to inner components including the core. At 
any pulse phase, low frequency emission comes from higher altitudes than high 
frequency emission. We have modeled the emission heights of pulse components 
of PSR B0329+54, and estimated field line curvature radii and particle Lorentz 
factors in the emission regions. 

Subject headings: stars: magnetic fields — radio continuum: stars — pulsars: 
individual (PSR B0329+54) — pulsars: general 



1. INTRODUCTION 

Pulsar radio emission is generally thought to be coherent radiation from the relativistic 
plasma streaming along the open magnetic field lines. The emission beam geometry has 
been widely attempted to interpret in terms of emission in purely dipolar magnetic field 
(e.g., Radhakrishnan & Cooke 1969; Goldreich & Julian 1969; Sturrock 1971; Ruderman & 
Sutherland 1975; Michel 1982; Blaskiewicz et al. 1991). 

Radhakrishnan and Cooke (1969) have proposed the rotating vector model (RVM) to 
explain the polarization-angle (PA) traverse (S-curve) of pulsar profiles. It has been fitted 
to average PA curves, and found to be quite successful (e.g., Lyne & Manchester 1988). The 
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three chief assumptions (Hibschman & Arons 2001) of RVM model are (i) the relativistic 
plasma flow is collimated by strong dipolar magnetic field lines, (ii) the observed radiation 
is beamed in the direction of tangents to the field lines, and (iii) the polarization angle is 
at a fixed angle to the curvature of field lines, as is the case of vacuum curvature radiation. 
Based on the assumptions (i) and (ii), we have attempted to locate the radio emission region 
in pulsar magnetosphere from which a distant observer receives the radiation. 

Statistical analysis of the distribution of pulse components within the pulse window 
indicates that the emission beam is conal (e.g., Mitra & Deshpande 1999; Kijak & Gil 2002). 
Rankin (1993) has indicated that inner pulse components are emitted at lower altitudes than 
the outer ones. Gangadhara and Gupta (2001, hereafter GG01) have estimated the emission 
height of different radio pulse components of PSR B0329+54, based on the aberration- 
retardation phase shift. They have also found that the inner components of pulse profiles 
originate from lower altitudes than the outer components. In the following paper Gupta and 
Gangadhara (2003) have then estimated the emission heights of 6 more pulsars using the 
same technique. Recently, the emission heights of all those pulsars have been re-estimated by 
Dyks, Rudak and Harding (2003, hereafter DRH03) by revising the aberration phase-shift 
relation given in GG01. Also, there are other claims for core and cone emission altitudes 
being different (e.g., Mitra & Rankin 2002) and roughly the same (e.g., Gil 1991). 

In this paper, we solve the viewing geometry for receiving the radio waves by a distant 
observer in § 2. By assuming the particle Lorentz factors, we show in § 2.2 that curvature 
radiation predicts the emission heights which are comparable to the radio emission heights 
in PSR B0329+54. In § 3, by utilizing the properties of curvature radiation and the emission 
altitudes estimated from aberration-retardation phase shift (DRH03), we estimate the field 
line curvature radii and particle Lorentz factors in the emission region. 

2. EMISSION-BEAM GEOMETRY 

Consider a magnetic dipole situated at the origin with magnetic dipole moment (m) 
aligned with the rotation axis (f2). In polar coordinates, the equation for a dipole-field line 
(Alfven & Falthammar 1963) is r = r c sin 2 9, where 9 is the magnetic colatitude and r the 
distance from the origin. The parameter r e is the field line constant, which is a distance from 
origin to the point of field line intersection with the equatorial plane (9 = n/2). In the case 
of an aligned dipole, r e = tlc for the last open field line, where tlc = Pc/2ir is the light 
cylinder radius, c the speed of light, and P the pulsar period. In the Cartesian coordinate 
system with z-axis parallel to fl, the position vector of an arbitrary point Q on a field line 
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is given by 

r c = r e (sin 3 9 cos0, sin 3 9 sin 0, sin 2 9 cos 9), (1) 
where <f> is the magnetic azimuth. 

Now consider the situation where the dipole (m) is inclined through an angle a with 
respect to fl, and rotated by phase 0' around z-axis, as shown in Figure 1. The position 
vector of the point Q on a magnetic field line, which is tilted and rotated, is given by 

r ct = A • r c , (2) 

where the transformation matrix A = R ■ I . The matrices for tilt (inclination) / and rotation 
R are given by 

cos a sin a \ / cos 0' — sin \ 

1 , R = I sin0' cos0' ] . (3) 
- sin a cos a / \ 1/ 

The matrix I produces clockwise rotation of the dipole around y-axis, and R counter- 
clockwise rotation around the z-axis. 

At Q, we find the tangent to the field line by evaluating b t = dr ct /d9, and the curvature 
by evaluating k t = dh t /ds , where bt = b t /|b t | is a unit tangent vector, ds is the arc length 
of the field line and \h t \ = r e ^/5 + 3 cos(2 9) sin 9j\f2. Therefore, the field line curvature 
radius is given by 

1 r c sin#[5 + 3cos(2#)] 3 / 2 
P ~W\~ 3V2 [3 + cos(2 9)} 

Since m is chosen to be parallel to z, the transformed magnetic dipole moment is given 

by 

m t = A • z = (sin a cos0', sin a sin0', cos a). (5) 

Hence the magnetic field of a dipole (Jackson 1975), which is inclined and rotated, is given 
by 

B t = 5 (^) 3 [3(f.m t )f-m t ], (6) 

where B is the surface magnetic field and tns ~ 10 km is the neutron star radius. It can 
be easily shown that B t satisfies V • B t = 0. 

The dominant magnetic field in the emission region is often shown to be consistent with 
being dipolar and to study the shape of pulsar radio beams (e.g., Narayan & Vivekanand 
1983; Lyne & Manchester 1988; Kramer et al. 1997). In a static dipole magnetic field 
approximation, the basic features of the pulsar magnetosphere can be understood. On the 
time scales of the order of pulse-phase bin, which is very small compared to the pulsar period, 
the rotating dipole may be treated, approximately, a static dipole. 
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2.1. Magnetic Colatitude and Azimuth of Emission Spot 

Consider the sight line n = (sin£, 0, cos( ), which lies in the xz-plane and makes an 
angle ( with respect to Cl, where ( = a + (3, and j3 is the angle of closest approach of sight 
line with respect to the magnetic axis. The half opening angle (T) of the emission beam is 
given by 

cos T = n ■ m t = cos a cos ( + sin a sin ( cos 0' . (7) 
If r is the angle between b t and m t then we have 

r ~ l + 3cos(2#) 
cosr = b t -m t = — , 8 
v/10 + 6 cos(2fl) 

In relativistic flow, radiation is beamed in the direction of field line tangent. So, at any 
instant, observable radiation comes from a spot in the magnetosphere where the tangent 
vector b t points in the direction of n. That is T m r at the emission spot. Therefore, using 
equations (7) and (8), we can find the relation for magnetic colatitude 9 as a function of T : 

I 

cos(2 9) = -(cos r Vs + cos 2 r - sin 2 r) , -7T < r < 7T . (9) 

o 

This equation is similar to the one given in GG01 (see, eq. 4) for 9. However, the equation (9) 
is superior compared to GGOl's equation as it represents a single physical solution when T 
changes sign. For T 1, equation (9) reduces to the well known approximate form 9 |r. 

Next, let k be the angle between n and b t , then we have 
cos k = n ■ b t 

= cos 2 T + (cos a sin ( cos 0'— sin a cos () sin T cos — sin ( sin 0' sin T sin 0. (10) 

Again to receive the radiation k ~ 0, therefore, we find the magnetic azimuth of the 
emission spot by solving n • b t = 1 and n x h t = : 

sin = — sin ( sin 0' esc T, 

cos = (cos a sin ( cos 0' — cos ( sin a) esc T. 



Therefore, we have 



sin C sin 6' 

= arctan — ^ ^— . (11) 

\cos(, sin a: — cos a sin( cos0'/ 



Note that is measured with respect to the meridional plane defined by Cl and m t . 
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2.2. Radio Emission Height 

Pulsar radio emission is generally believed to be coherent curvature radiation by sec- 
ondary pair plasma streaming along the dipolar magnetic field lines. The characteristic 
frequency of curvature radiation, at which the emission peaks, is given by (Ruderman & 
Sutherland 1975): 

3c 7 3 
Air p 

where 7 and p are the Lorentz factor and the radius of curvature of trajectory of a rela- 
tivistic particle, respectively. Since particles closely follow the curved dipolar field lines, the 
curvature of particle trajectory can be approximated with the curvature of field lines. 

For the given v and 7, equation (12) predicts the value of p required. Then using equa- 
tion (4) we can find the field line constant r c with the help of 6 em obtained from equation (9). 
The azimuth cm of the emission spot can be determined by using equation (11). Hence, 
by having defined (r e , 6 em , em ), we can estimate the emission height r em of radio waves of 
frequency v with respect to pulse phase 0' using equation (2). 

Based on the aberration-retardation phase-shift of pulse components in observed profiles, 
Gangadhara and Gupta (2001) have estimated the emission heights of pulse components in 
PSR B0329+54 at frequencies 325 MHz and 606 MHz. Recently, Dyks, Rudak and Harding 
(2003) have revised the aberration phase-shift relation given in GG01, and re-estimated the 
emission heights. In Figure 2, we have plotted the revised emission heights (marked points 
with error bars). 

For modeling the emission heights in PSR B0329+54, we adopt (a, (3) = (30°, 2.1°) given 
by Rankin (1993). We find that the curvature emission predicts the emission heights (solid 
and dashed line curves in Figure 2) which are comparable to the observed ones (DRH03) 
if the secondary particle Lorentz factor 7 ps 340 and 390 for the emissions at 325 MHz 
and 606 MHz, respectively. Figure 2 shows that at any given frequency the emission near 
the pulse center (core) comes from lower heights compared to outer conal components, in 
agreement with the results derived from observations (marked points). At any pulse phase, 
low frequency emission comes from higher altitude than the high frequency emission, in 
agreement with the well known radius-to-frequency mapping (RFM) (e.g., Cordes 1978; 
Phillips 1992). We find RFM mapping is not uniform across the pulse window, i.e., it 
is more pronounced in the case of outer cons (say, cones 3 & 4) than in the inner cones 
(say, cones 1 & 2). Mitra and Rankin (2002) have also made a similar proposition from the 
study of 10 selected pulsars. The emission height difference 5r em between the two frequencies 
increase progressively with respect to pulse phase, and the marked points in Figure 2 indicate 
5r em ps 14 + 0.046 0' 2 + 0.003 0' 4 km, where the pulse phase 0' is in degrees. 
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2.3. The Polar Cap 

The polar cap boundary is defined by the last open field lines for which k t • z = at the 
light cylinder. Therefore, the magnetic colatitude 9\ { = 9, at which k t ■ z vanishes, is given 
by 

, , -3ai 2 -a 2 V8ai 2 + a 2 2 

COS(2 lof ) = q 2 i 2 ' -7T < < 7T 13 

9 a\ z + a2 2 2 

where a\ = sin a cos0 and a 2 = 3 cos a. Next, for the range of n/2 < <p < Stc/2, the 
colatitude is given by it — 9\ ol . The angle (rj) between f ct and z is given by 

cosrj = z • r ct = cos a cos 9 — sin a sin 9 cos0. (14) 

If Vm = V at the light cylinder for last open field line then we have |r ct | sin^f = r LC , and 
the last open field line constant 

r e ,iof = r LC esc 2 6>iof esc r)i of . (15) 
Hence the radial position of foot of a last open field line from the magnetic axis is given by 

s P = r NS 9 P , (16) 

where 9 V = arcsin (y/r^sj r e,iof ) is the colatitude of foot of last open field line. 

To plot polar cap consider the Cartesian coordinate system: (xb, jb, zb) such that the 
axis zb is parallel to m t and xb lies in the meridional plane. The coordinates of foot of a 
last open field line is given by 

(xb ,Vb) = (Spcos0,s p sin0). (17) 

Using a = 30°, the polar cap of PSR B0329+54 is plotted in Figure 3. It is nearly elliptical 
with the dimension of 164 m in the xb direction and 171 m in the direction of yB, in agreement 
with the polar cap shapes proposed by Biggs (1990). 

For a given pulse phase 0', using equations (9) and (11), we can find 9 em and cm of 
the emission spot. The radial location of foot of field lines, from which a distant observer 
receives the radiation, is given by 

(18) 

Thus, by knowing s cm and </> em one can specify (xb, Ub) of the foot of field lines, which 
direct the radiation beam towards the observer. The values of s em estimated from the 
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conal aberration-retardation phase-shift (marked points with error bars in s em ) are plotted 
in Figure 3. The model (dashed and solid line) curves represent the locus of foot of those 
field lines, which are accessible to a distant observer in 325 MHz and 606 MHz observations, 
respectively. We find the low frequency emission is received from the field lines which lie 
closer to the magnetic axis than the high frequency ones. The curve (solid line) indicates 
increasing conal ring spacing (5s em ) between the successive rings from inner (cone 1) to outer 
(cone 4) on polar cap: ~ 14 m between the 1 st ring and the 2 st ring, 16 m between 2 nd and 
3 rd , and 19 m between 3 rd and 4 th . The dashed line curve for 325 MHz emission also indicate 
the similar ring spacings, except they are less by ~ 0.5 m. These results tend to support 
the model of concentric rings of sparks produced in the vacuum gap region just above the 
neutron surface (e.g., Ruderman & Sutherland 1975; Gil & Sendyk 2000). 



3. DISCUSSION 

In the previous section, to explain the emission height of components in PSR B0329+54, 
we selected the Lorentz factor 7 of about 340 and 390 for the emissions at 325 MHz and 
606 MHz, respectively. It means p and 7 assume constant values across the pulse window. 
However, it may not be reality, as the model heights significantly depart from the observed 
ones (marked points) as indicated by Figure 2. On other hand, one can turn around the 
calculation by accepting the observed emission heights ask for the required 7 and p. In 
column 2 of Tables 1 & 2, we have given the phase location of conal components in the 
absence of aberration-retardation phase shift (see, eq. 11, GG01), i.e., in corotating frame. 
Next, the columns 3 and 4 give the magnetic colatitude 9 and the longitude of emission 
spot, estimated using equations (9) and (11). The small values of 9 < 6° implies the conal 
emission mostly comes from the vicinity of magnetic axis. By having known r em and 9 em , we 
estimated the radius of curvature p using equation (4), and given in column 5. It shows inner 
cones are emitted from the region of smaller curvature than outer ones. Next, for a given 
frequency v and known p, equation (12) allows to estimate the 7 of secondary particles, as 
given in column 6. It shows particles with slightly higher 7 move on the field lines associated 
with outer cones than inner ones. 

The emission beam by relativistic particles is centered on their velocity, which is roughly 
parallel to the field line tangent b t , and has a radial angular width of I/7. Consequently, 
observer tends to receive the radiation even when n does not perfectly align with b t , which 
can lead to a spread in the emission height as well as in pulse phase of component peak 
location. So, we estimated the spread in emission height of conal components by allowing 
T in equation (9) to vary by I/7. However, we find the spread is minimal and lies within 
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the error bars. For example, at 325 MHz, we find 9 km for the 1 st cone and 15 km for 
the 4 th cone. At 606 MHz, it is 6 km for the 1 st cone and 10 km for the 4 th . For other 
cones the spread lies in between these values. We also estimated the spread in phase of 
component peaks, and found to be less than 0.4° for all the components. If we compare the 
emission heights (DRH03) along with these spreads, we find the emission regions of conal 
components at the two frequencies are well separated: it is 14 km in the case of 1 st and 214 
km for the 4 th . Further, the pulse phase locations of conal components are different at the 
two frequencies, as indicated by column 2 in Tables 1 & 2. So, the observer tends to see the 
components separated in phase as well. Note that the received radiation is maximum only 
when n perfectly aligns with b t . Hence we may conclude that the observer tends to receive 
the radiation at the two frequencies from two distinctly located regions in the magnetosphere. 

Using these values of 7 and p, we estimated the power radiated by an electron or a 
positron by curvature emission (see, eq. 25a, Ruderman & Sutherland 1975), and given in 
column 7. It shows the single particle emission is maximum in the region of inner cone 
than in the outer cone regions. However, we know from the estimates made by Ruderman 
and Sutherland (1975) that the incoherent superposition of single particle radiation cannot 
explain the pulsar radio emission. The very high brightness temperature (10 25 -10 30 K) of 
pulsars lead to the conclusion that it must be coherent. The coherence due to the bunching of 
plasma particles in the magnetosphere was proposed much earlier (e.g., Pacini & Rees 1970; 
Sturrock 1971). The coherence due to bunching goes as inverse of wavelength of radio waves, 
therefore, the emission by bunches at low frequency tends to dominate the high frequency 
emission. 

Thus, we have developed a method for estimating the possible values of 7 and p needed 
for the emission of radio waves at a given frequency. By matching the coherent curvature 
emission with the observed pulsar fluxes, one may be able to estimate the plasma density 
and the coherency (bunching) factor exist in the radio emission region of pulsars. For v > v v 
strong coherent radiation occurs, where v v is the plasma frequency (Ruderman & Sutherland 
1975). 

While estimating the aberration-retardation phase shift of conal components, the core 
(reference) phase was set to zero (GG01; DRH03). However, it is possible that the core is 
also emitted at a finite height from the neutron star surface, as indicated by Figure 2. We 
fitted a 4 th de gree polynomial to the emission heights (marked points in Fig. 2) of 8 conal 
components. At zero phase, the fit indicated emission height of 26 ± 50 km for the 325 MHz 
points and 12 ± 45 km for the 606 MHz points. These heights may be treated as minimum 
values for the core emission height, as the component heights themselves have been estimated 
by setting the core phase to zero. Though it is important to consider the emission height 
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of core, we believe our results will not change significantly by the core height, as aberration 
and retardation are minimal at that height. 

4. CONCLUSION 

By assuming the curvature emission by relativistic sources is beamed in the direction 
of tangent to the dipolar magnetic field lines, we have resolved the pulsar radio emission 
geometry by solving the viewing geometry in an inclined and rotating dipolar magnetic 
field. Due to the geometric restrictions, a distant observer tends to receive outer conal 
components from higher altitudes than inner ones including the core. Further, we find low 
frequency emission comes from higher altitudes than high frequency emission, in agreement 
with the radius-to-frequency mapping (RFM). We find RFM is more pronounced in the 
outer conal components than in the inner ones. Using the emission heights obtained from 
aberration-retardation phase shift, we have estimated the probable values of Lorentz factor 
and field line curvature radii in the emission region. 

I would like to thank Don Melrose and Simon Johnston for illuminating discussions, 
and J. M. Rankin, J. A. Gil, Y. Gupta, V. Krishan, and C. Sivaram for comments. I thank 
anonymous referee for useful comments. 
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T able 1: Parameters related to radio emission from PSR B0329+54 at 325 MHz 

Cone ^ 8 4> p/r LC 7 P/Pi f 

No. (deg) (deg) (deg) 

1 5.0 ± 0.26 2.22 ± 0.07 -53.0 ± 1.6 0.15 ± 0.00 286 ± 3 1.00 ± 0.08 

2 8.4 ± 0.19 3.21 ± 0.06 -67.8 ± 0.6 0.23 ± 0.00 329 ± 2 0.76 ± 0.03 

3 12.0 ± 0.27 4.35 ± 0.09 -76.5 ± 0.5 0.31 ± 0.01 363 ± 2 0.62 ± 0.03 

4 17.0 ± 0.80 5.99 ± 0.26 -83.9 ± 1.0 0.33 ± 0.01 370 ± 5 0.60 ± 0.06 



t Pi = 1.16 x 1(T 16 erg s" 1 



T able 2: Parameters related to radio emission from PSR B0329+54 at 606 MHz 

Cone ft 4> p/r LC 7 P/Pj 

No. (deg) (deg) (deg) 

1 4.9 ± 0.25 2.19 ± 0.07 -52.4 ± 1.6 0.12 ± 0.00 328 ± 3 1.00 ± 0.07 

2 7.8 ± 0.18 3.02 ± 0.06 -65.8 ± 0.6 0.21 ± 0.00 391 ± 2 0.70 ± 0.03 

3 11.0 ± 0.45 4.02 ± 0.14 -74.5 ± 1.0 0.26 ± 0.01 419 ± 5 0.61 ± 0.05 

4 15.7 ± 0.59 5.56 ± 0.19 -82.3 ± 0.8 0.26 ± 0.01 420 ± 5 0.61 ± 0.05 



t Pi = 3.06 x 1(T 16 erg s" 1 
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Fig. 1. — Coordinate system to describe the inclined and rotating magnetic dipole. The 
rotation phase 0' is measured from the fiducial plane (x-z) in the counter-clockwise direction 
around Cl. The magnetic colatitude 9 is measured from m t , and the magnetic azimuth (f> 
measured counter-clockwise around m t from the meridional plane defined by Cl and m t . 
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Fig. 2. — Radio wave emission heights in PSR B0329+54 with respect to rotation phase 0': 
solid and dashed line curves are for the emissions at 606 MHz and 325 MHz, respectively. The 
emission heights estimated from aberration-retardation phase shift (DRH03) are superposed 
for comparison: the points marked with o and • are for the components at 325 MHz and 
606 MHz, respectively. 
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Fig. 3. — Polar cap of PSR B0329+54. Solid and dashed line curves represent the locus of foot 
of field lines associated with the emissions at 606 MHz and 325 MHz, respectively. The values 
of s em estimated from aberration-retardation phase shift are superposed for comparison: the 
points marked with o and • are for 325 MHz and 606 MHz, respectively. 



